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Abstract 

The ECO method and the theory of Catalan-like numbers intro- 
duced by Aigner seems two completely unrelated combinatorial set- 
tings. In this work we try to establish a bridge between them, aiming 
at starting a (hopefully) fruitful study on their interactions. We show 
that, in a linear algebra context (more precisely using infinite matri- 
ces), a succession rule can be translated into a (generalized) Aigner 
matrix by means of a suitable change of basis in the vector space of 
one-variable polynomials. We provide some examples to illustrate this 
fact and apply it to the study of two particular classes of succession 
rules. 

Keywords - Succession rules, Catalan-like numbers, Pascal matrix. 

1 Introduction 

The ECO method was founded in the 90's by a group of researchers, 
including Pinzani, Barcucci, Del Lungo and Pergola [BPLPP1 iBPT/PPIj . 
It consists of a purely combinatorial way of constructing the objects of a 
given class in such a way that, if the construction is sufficiently regular 
and recursive, enumeration follows by more or less standard methods of 
combinatorial analysis. More precisely, one starts by partitioning a class of 
objects according to their size (suitably denned). The goal is then to perform 
a sort of local expansion on each object of a given size, thus producing all 
the objects of the successive size exactly once. Therefore a single object 
produces a set of new objects according to some parameter. Typically, if 
such a construction is regular enough, one can encode it using a succession 
rule |W11 IW2j . which is a purely formal system, generally expressed as 
follows: 

(a) (1) 
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Here the letters denote positive integers, (a) is called the axiom and 
(k) (ei(k)) ■ ■ ■ (efc(fc)) is the production of (k). A succession rule can 
be represented by means of its generating tree, which is, by definition, 
the infinite, rooted, labelled tree whose root is labelled a (like the axiom) 
and such that every node labelled k produces k sons, labelled respectively 
ei(fc), . . . , efc(fc). One of the main enumerative information provided by a 
succession rule is the numerical sequence of the cardinalities of the levels of 
the generating tree associated with the rule: we will refer to such a sequence 
as the numerical sequence determined by the succession rule. 

The basic reference for the ECO method is BDLPP , in which many 
examples can also be found. 

The importance of succession rules as a tool for the ECO method has 
lead to several investigations to get a better mathematical insight on them. 
In |FPj the authors define the concept of rule operator, thus translating a 
succession rule into a linear operator on one-variable polynomials. In DFR 
every succession rule is associated with at least two infinite matrices: the 
production matrix, which is essentially the matrix of the related rule operator 
with respect to the canonical basis of polynomials (s") n gN, and the ECO 
matrix, whose (n, /c)-entry is, by definition, the number of nodes labelled k 
at level n in the corresponding generating tree. We point out that a few 
instances of the notion of production matrix appeared for the first time in 
|Wlj under the name of "transfer matrices" . 

Another combinatorial theory dealing with infinite matrices is Aigner's 
theory of Catalan-like numbers [MJ EH |Zj. The basic idea of [H] is 
to characterize those sequences for which the determinants of the Hankel 
matrices of order are all equal to 1. It is shown that such sequences appear 
as the first column of certain infinite matrices, called admissible matrices. 
These numbers are referred to with the name of Catalan-like numbers. The 
reason for this name lies in the fact that Catalan numbers are the unique 
sequence whose Hankel determinants of orders and 1 equals 1. 

In |A21 IA3j Aigner extends this theory by considering a more general 
kind of matrices, which we will rename Aigner matrices (instead of the infe- 
licitous name "recursive matrices" used in |A3j ) . Generalizing the previous 
definition, we will call Catalan-like numbers every sequence appearing as the 
first column of an Aigner matrix. 

The aim of our work is to provide a "vocabulary" to translate the ECO 
method into Aigner's theory, and vice versa. Such a vocabulary turns out to 
be based on linear algebra tools, consisting of a suitable change of basis in 
the vector space of one-variable polynomials. What we hope to show in this 
paper is that the two theories under consideration are, in some algebraic 
sense, the two sides of the same medal, which is quite surprising if we think 
of the very different starting points, and combinatorial meanings, of such 
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theories. 

After a brief survey of the notions we need from the ECO method and 
Aigner's theory, we provide the main linear algebra tools to be used in the 
sequel. In particular, we define the Aigner basis in the vector space of one- 
variable polynomials and prove some of its properties. Next we introduce 
what we call the fundamental change of basis, which is the key ingredient 
to accomplish our project, and propose some examples to illustrate our ap- 
proach. The final part of the work is devoted to the study of two particular 
cases, for which we are able to fully describe how to switch from one theory 
to the other. In the last section we give some hints for possible, future works. 

Before starting, a few words concerning our notations. We have chosen 
to use N to denote natural numbers (zero included), whereas N* is the 
set of natural numbers without zero. The symbol x is used to denote the 
operator of multiplication by x, in order to distinguish it from the symbol 
x, used as a polynomial variable. The symbol T denotes the transpose of a 
matrix. The last remark concerns the way we have chosen to index the lines 
of our matrices. Classically, the lines of an Aigner matrix are indexed by N, 
whereas, in an ECO matrix, the first column is usually column 1 (so that 
columns are indexed by N*). There are clear combinatorial and algebraic 
reasons for this: the n-th column of an ECO matrix gives the distribution 
of label (re) in the generating tree, whereas the scalar product of the n-th 
and the m-th rows of an admissible matrix gives the first element of its 
(n + m)-th row. Unfortunately, keeping both these conventions would result 
in a purely formal, but greatly inelegant, variation of our results: namely, 
the fundamental change of basis (which is degree-preserving in our theory) 
would translate x n into a polynomial of degree re — 1 (which is p n ~i(x), 
according to the notations of section |HJ). To avoid this we have preferred to 
use N* as the set of indices for the lines of Aigner matrices. In this way 
the (nice) characteristic algebraic property of admissible matrices becomes 
a little bit difficult to read, but our theory can be described much more 
elegantly. 

2 Preliminaries on ECO and Aigner matrices 

In this section we report the main facts concerning the two combinatorial 
theories we are going to compare. 

Consider a succession rule as in (pQ) . Instead of representing it by means 
of a generating tree, one can choose linear algebra notations. In the vector 
space of one-variable polynomials, define the linear operator L = Lq on the 
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canonical basis (x™) ng N as follows: 

L(l) = x a , 
L(x k ) = x ei{k) + . . . + x ek{k) , 
L(x h ) = hx h , if (h) is not a label of fl 

L is called the rule operator associated with Q (see |FPj ): it bears all 
the enumerative properties of a succession rule and allows to express such 
properties using algebraic notations. For example, if (/ n )neN is the numerical 
sequence determined by f2, then we can find f n using the rule operator L of 
f2 as follows: 

/„ = [L n+ \l)] x=l . 

Throughout the present work, we will always deal with a special case, 
namely we assume that deg L n (l) = n. From a combinatorial point of view, 
this means that the set of the labels of 0, is N* and the maximum label 
among those produced by (k) is (k + 1). 

The infinite matrix P = Pq representing L with respect to the canon- 
ical basis (x n ) n gN is called the production matrix of f2. Such matrices are 
extensively studied in |DFR| : here we only recall some of their properties. 

Let Ap be the infinite matrix whose n-th row vector is given by n T P n ~ 1 
(where u T = (1 ... . . .)). Then Ap describes the statistic given by the 
distribution of the labels at the various levels of the generating tree related 
to Q. Namely, the (n, k) entry of is the number of nodes labelled k at level 
n of the generating tree of f2. Ap is called the ECO matrix associated with P 
(or with £1), and is also characterized by the matrix equality DAp = ApP, 

...\ 

... 

1 ... 
1 ... • 
... 

V : ■■) 

In terms of the production matrix P, the sequence (/ n )neN is nothing else 
than the sequence of the row sums of the associated ECO matrix. The 
ordinary and exponential generating functions of Q are given, respectively, 
by fp(t) = u T (I — tP)~ 1 e and Fp(t) = u T exp (tP)e, where e is the column 
vector (111 ...1 •••) T and exp (X) denotes the usual matrix exponential 
of the (infinite) matrix X. In section we deal with some examples of the 
theory we are going to develop; in describing such examples we also consider 
the production and ECO matrices of some classical succession rules. The 
reader is invited to have a look to those examples in order to be introduced 
to these concepts. 
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In a recent series of nice and well- written papers |A1| IA2| IA3j , Martin 
Aigner has developed a new theory to deal with numerical sequences some- 
how linked to the sequence of Catalan numbers. One of the main tool of this 
theory is a particular class of infinite (triangular) matrices, called admissible 
matrices in |Alj and renamed with the infelicitous term recursive matrices 
after |A3j . In the sequel we will use the following terminology. 

Consider a lower triangular matrix A = (a n fc)n,fceN* with main diagonal 
equal to 1. A is called an admissible matrix whenever, for every n, m € N*, 
the (ordinary) scalar product of the n-th and the m-th rows of A gives the 
first element of the (n + m — l)-st row; in symbols: 

^ ^ O'nkO'mk = ^(n+m— 1)1 • 
fc>l 

More generally, we define A to be an Aigner matrix when there exists 
a sequence of nonnegative integers (T n ) nS N* such that T\ = 1 and T n \T n+ \ 
for which we have 

^] O'nk&mkXk = a (n+m-l)l- (2) 
fe>l 

Obviously admissible matrices correspond to Aigner matrices for which 
T n = 1. 

If A is an Aigner matrix, the sequence (a n i) nS N* is called the sequence 
of Catalan-like numbers associated with A. For an extensive study of the 
algebraic and enumerative properties of Aigner matrices the reader is re- 
ferred to |A1[ IA3j and to the further items cited in the references of the two 
papers. Here we recall only those results which we need for our purposes. 

Proposition 2.1 ( X4i[l^ffi ) An Aianer matrix A = (a n fc)n,fceN* is uniquely 
determined by the two sequences (a(n+i)n)neN* an d (^n)neN* • Conversely, to 
every pair of sequences (6 n )neN* and (T n ) ne N* of real numbers there exists 
an (and therefore precisely one) Aigner matrix A = (a n k)n,k£N* associated 
with (T n ) nG N* and such that a( n +i) n = b n for all n. 

Proposition 2.2 (I A 1\ \A^ ) Let A = (a n k) n ,keN* be an Aigner matrix as- 
sociated with (T n ) neN ». Set si = a%\, s n = a(„+i)„ - a n (n-i) an d tn = 
for n > 2. Then we have 

a n = 1, aik = (k > 1) 

«nfc = d(„_l)(fc_l) + Sfca(„_!) fe + tfc+lO(n-l)(Jfe+l) ( n - 2 )- ( 3 ) 

Conversely, if a n k is given by the recursion |pj), then A = (a n fc) n ,fceN* 
is an Aigner matrix with T n = ti ■ . . . ■ t n and ar n +x\ n = s± + . . . + s n . 
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Setting a = (s ra ) n eN* ; t = (t n ) n >2, we say that A = A a ' T is the Aigner 
matrix of type (a, r) when recursion (j3J) holds for its entries. 

It is possible to write (j2j) in a compact matrix form. Define the diagonal 
matrix T and the infinite Hankel matrix of the sequence (a n i) ne N* as follows: 

/ an a 2 i 031 • • • \ 

021 031 a 41 ' • • 
«31 041 a 51 ' - - 

Then (J2J) can be written as ATA T = H. More precisely, we have the 
following characterization. 

Proposition 2.3 (]Aty ) A is an Aigner matrix if and only if ATA T = H 
with T n 7^ for all n > 2, T\ = 1. The sequences a and r are then given as 
in proposition \2.'lA 

Moreover, if we denote by H n the n-th Hankel matrix of a sequence 
(a n ) n£ N* (which is, by definition, the submatrix of H consisting of rows and 
columns 1 to n), we have the following corollary. 

Corollary 2.1 ( ]A Sf ) A sequence (a n ) ng N* is Catalan-like if and only if 
\H n \ / for all n > 1. 

3 Linear algebra tools 

In this section we give some (elementary) linear algebra tools which will 
be necessary in the sequel to settle the stated analogy between ECO method 
and Aigner's theory of Catalan-like numbers. 

In the vector space of one- variable polynomials over the real field (to be 
denoted R[x]) we define the following polynomial sequence: 

Po(x) = 1, 

Pn (x) = x{x-l) n ~ l , Vn>l. 

It is clear that deg p n (x) = n, so that (j? n (^))neN constitutes a basis for 
the vector space R[x]. We call such a basis the Aigner basis of R[x]. 

Remark. We recall that the polynomial p n (x) has a very important com- 
binatorial meaning: it is the chromatic polynomial of a tree having n vertices 
(see, for example, B ). However, this fact will not be used in this paper. 

It is well known (see, for instance, |A0| ) that, for any basis of R[x], there 
exists a unique differential 1 linear operator mapping the n-th element of the 

i.e., it maps a polynomial of degree n into a polynomial of degree n — 1. 
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basis into the (n — l)-st one. According to |FPj . we call factorial derivative 
operator the linear operator: 

T : R[x] — > R[x] 
■ Po(x) — > 0, 

: Pn(x) >p n -i(x), n > 1. 

The factorial derivative operator can also be nicely performed on the 
canonical basis of R[x]. 

Proposition 3.1 For any n G N, T(x n ) = 1 + x + . . . + x n_1 = YJk=l x k ■ 

Proof. For the first values of n, we have T(l) = 0, T{x) = 1, T(x 2 ) = 
T(pi(x) +P2(a^)) = 1 + x. By induction, suppose that T{x n ) = 1 + . . . x^ 1 . 
Observe that we have the trivial equality T((x — l)p n (x)) = T{p n+ \(x)) = 
Pn(x), whence T{xp n {x)) = T(p n (x)) + p n (x). By linearity we then have: 

T(xp(x)) = T(p(x)) +p(x). 

Therefore, in the case p(x) = x n , we get: 

T{x n+1 ) = T{x-x n ) = T{x n ) + x n 
= l + .-. + x^ + x", 

which is the thesis. ■ 



Some properties of the factorial derivative operator are recorded in |FPj . 

The Aigner basis has a nice behavior with respect to the usual multi- 
plication operation on polynomials. 

Proposition 3.2 Let n,m>l. 

1. p n (x) ■Pm(x) = Xp n+m -i(x); 

2. x k p n (x) = J2h=o [h)Pn+h(x); in particular, xp n (x) = p n +l(x) +Pn{x); 

3. denoting by x _1 the linear operator defined by x _1 (p(x)) = p ^~ p ^ 
(so that x -1 is the usual difference quotient operator), it is x~ fc p n (x) = 
Sft=o( — l) n ~ k ~ h ( n Z^\ fl )Ph(x); in particular, setting k = 1, we have 

X _1 p n (x) = ^=d(-l) n_/l_1 Pft.(a;) = p n -l{x)-Pn-2{x)+Pn-?,{x)- ■ ■ ■■ 

Proof. 

1. p n (x) ■ p m {x) = x(x - l) n ~ l ■ x(x - l)™" 1 = x- x(x - l) n + m ~ 2 = 
xp n+m ^i(x). 
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2. We have immediately 



p n+ l(x) +p n (x) = x(x -l) n X (x - 1 + 1) = Xp n (x). 

Then, by induction, we get 

X k+1 p n (x) = X ■ X k p n (x) = X^ \ V\Pn+h{x) 

h=0 ^ ' 

= E ( tj^+h+lix) +p n +h{x)) 



h=l 

fe+1 



E{ k+ h %+^)- ( 4 ) 



/i=0 



3. Just recall that x k is the inverse of the operator x fc . Then the thesis 
is obtained by simply inverting the combinatorial sum in @. ■ 

We close this section by stating a technical result, useful in the com- 
putation of the powers of the factorial derivative operator T, which can be 
proved by induction. 

Proposition 3.3 For n, k G N, we have: 

h=0 v J 
Corollary 3.1 1. [T k {x n )] x=x = (™); 
2. [T k (x n )] x=0 = ( n k Zl). 

4 The fundamental change of basis 

For our purposes, we slightly modify the definition of ECO matrix given 
in |DFE,j . namely we suppose that the n-th row of F (with n G N*) gives the 
distribution of the various labels at level n — 1 of the generating tree of the 
rule; thus the (n, k) entry of F is the number of nodes labelled k at level n—1. 
Let Q be a succession rule as in (^) and suppose that F = (/ n fc)n,fcgN* is the 
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ECO matrix associated with Q, as explained above. We denote by (5 n {x) the 
polynomial canonically associated with the n-th row of F, namely: 

n 

Pn(x)=Y,fnkX k (n>l). (5) 

k=l 

As we have already remarked, we will always assume that deg (3 n {x) = n. 
Now expand the polynomials f3 n (x) in terms of the Aigner basis, thus 
obtaining 

n 

0n{x) = ^2a nk p k (x). (6) 

k=l 

Clearly, the coefficients f n k and a n k are intimately related. In particular, 
the following, simple result shows that this setting is the right one to achieve 
our project. 

Proposition 4.1 The succession (r„) ne N* = (Sjt=i 7nfc)neN* of the row 
sums of F is equal to the succession (a n i) nG N* given by the first column of 
the matrix A = (o n fc)n,fceN* • I n symbols: 

n 

Q"nl = ^ ~] fnk- 
k=l 

Proof. It is obvious that r n = f3 n (l), whence YHt=i a nkPkiX) = r n . By 
the definition of the Aigner basis, it is Pfc(l) = 0, for k > 1, and = 1, 

and so r n = a n i, as desired. ■ 

The results obtained so far can be naturally expressed also in matrix 
notation. Specifically, it turns out that the fundamental change of basis 
described in Q is represented as the multiplication on the right by the 
Pascal matrix. In other words, if P = ((^)) n fcgN is the usual Pascal matrix, 
and F = (/ n fc)nfeeN* a given ECO matrix, then the associated matrix A = 
(«nfc)n,fceN* can be expressed as follows: 

A = FP. 

So the Pascal matrix P is the matrix of the change of basis from (x n ) nS N 
to the Aigner basis ip n (x)) n e'N- 

Now cards are laid on the table: changing the canonical basis into the 
Aigner basis is the "linear algebra" way to switch from the ECO method to 
Aigner's theory. At this stage, the following, very natural question can be 
asked: 
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1) for which ECO matrices F does it happen that the matrix A is an 
Aigner matrix? 

Regarding things the other way round, one can start with an Aigner ma- 
trix A and perform the inverse change of basis (from p n (x) to x n ). Therefore 
the previous question can be inverted: 

2) for which Aigner matrices A does it happen that the matrix F is an 
ECO matrix? 

These two problems seems to be rather difficult to be tackled in their full 
generality. In the rest of the paper we will mainly focus on special examples 
to hopefully illustrate the interest of our approach. Finally, we will consider 
two particular classes of ECO matrices (namely, those arising from the so- 
called factorial succession rules and differential succession rules), giving for 
them a complete answer to question 1. 

Before closing this section, we record some further notations and results 
which will be useful in the sequel. 

If L is the rule operator associated with f2, then for the polynomials 
Pn(x) in © we clearly have 



Applying L means to shift from row n to row n + 1 in the ECO matrix 
associated with O, whence: 



The coefficients and a n k in (J5J) and © can be expressed in linear 
algebraic terms, as the following proposition clarifies. 



f3 n (x) = L n (l). 



L{[3 n {x)) =/3 n+1 (x). 



(7) 



Proposition 4.2 



1- fnk 



f (&(*)) 



J x=0 



2. a n k 



'))] 



2 = 0' 



The easy proof is left to the reader. 



Corollary 4.1 For any n, a (n+1)n = /( n+ i)„ + n/ (n+1)(n+1) . 



Proof. From the above proposition we have 



a (n+l)n 



[T n (p n+1 (x))] x=0 



J2f(n + l)k[T n (x k )} 



x=0- 
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Now, recalling corollary I3.1[ we get immediately: 



a, 




k>l 



n+l)(n+l)> 



as desired. ■ 

To conclude, we prove the following proposition, concerning the behavior 
of a rule operator when applied to the Aigner basis. 

Proposition 4.3 For any rule operator L and for any n > 2, we have: 



(This is due to the fact that [L(x k )] x= i = k = [D(x k )] x= i). Then we get 



[L(p n (x))] x =i = [D( Pn (x))] x=l = [Dx(x - l)"" 1 ]^! = 0. ■ 

5 Detailed examples 



In this section we will provide a detailed analysis of the effects of the 
fundamental change of basis in the case of a well-known succession rule 
determining Catalan numbers. Then some other examples will be dealt with; 
for them, we will only state the main facts (without giving proofs), however 
our results can be checked out by a direct computation or by applying the 
theory we are going to develop in the final part of our work. 

5.1 Catalan numbers 

Let us consider the succession rule 



[L{Pn{x))] x =l = 0. 



Proof. It is easy to see |FPj that, for any rule operator L, it is 



[L(p(x))] x=1 = [D(p(x))] x=1 . 



immediately: 



Q:t (l)-(2) 

( (k) - (2) (3) (4) ...(*) (* + 1) 
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defining Catalan numbers 1,1,2,5,14,42,132,. . . (see, for example, lBDLPP]). 
The first lines of the ECO matrix associated with 0, looks as follows: 



F 
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1 
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" ^ 
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1 
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1 
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1 
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14 


14 


9 


4 


1 
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42 


42 


28 


14 


5 


1 


• 









132 


132 


90 


48 


20 


6 


1 • 




V 





















Applying the fundamental change of basis, one gets the following matrix: 

\ 



A 



(I 

110 

2 3 10 

5 9 5 1 

14 28 20 7 1 

V: ! : i : 



/ 



Is A the (unique) admissible matrix for Catalan numbers found in |Alj ? 
The first fact suggesting that it could be so follows from the application of 
corollary 14.11 It is well known |FPj that the entries of F are the so-called 
ballot numbers, namely: 



nk 



n 



1 f2n - k - 2 

n — k 



1 



(n, k > 1). 



Therefore it follows immediately that, in A, we have: 

a {n+l)n = /(n+l)n + n f(n+l)(n+l) 



n — 1 (n\ n fn — 1 

n VI/ +n 'n{ 



n — 1 + n = 2n — 1, 



which agrees with the (n + 1, n) entry of the admissible matrix for Catalan 
numbers. Obviously, this is not enough to conclude, however it is in fact a 
strong hint. To get to the desired result we need to show, for example, that 
the coefficients a n ^ obey the following recursion (deduced from |Alj ): 

a (n+l)k = a n(k~l) + ^ a nk + d n (k+l)- 

Expressing the a n fc's as in proposition 14.21 we then find: 

[T k (p n+1 (x))] x=0 = [(T k ~ l + 2T k + T k+1 )(P n (x))] x= o, 
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whence, recalling corollary 13.11 



eQ_J)/(^=e(Q_2) +2 

ft=i v 7 /i=i v v 7 



t |]/* (8) 



The sum of binomial coefficients in the r.h.s. of (jHJ) can be easily sim- 
plified (using well known properties of the Pascal matrix) to obtain: 

E u_ = E rt j^- ( 9 ) 

To prove equality Q we make use of the structural properties of the 
ECO matrix F, namely the recursion: 

f(n+l)h = fn(h-l) + fnh + • • • + fnn 

n 

= ^ 1 fni ■ 



i=h-l 



Replacing in the l.h.s. of © and interchanging the order of the summa- 
tions when necessary, we get: 



E^T)^ = E(fc_j) ( E fni) 
h=l v 7 h=l v 7 \i=h-l / 

- e(e(::;))/-ec:>- 

i=0 \h=l v 7 / i=0 v 7 

which is an identity (/ n o = by convention). Therefore, we have formally 
proved that switching from the canonical basis to the Aigner basis translates 
the ECO matrix F of Catalan numbers into the (unique) admissible matrix 
A of Catalan numbers. 

5.2 Motzkin numbers 

We can use the same approach to deal with Motzkin numbers. Consider 
the succession rule 

f (1) 

n : < (1) -w (2) 

k (fc)^(l)(2)(3)---(fc-l)(fc + l) 
Q determines the sequence of Motzkin numbers ( (BDLPP] ). and its ECO 
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matrix is the following: 



F 



( 1 




















1 














1 





1 











1 


2 





1 








3 


2 


3 





1 





6 


7 


3 


4 





1 



V 



/ 



Applying the fundamental change of basis (or, equivalently, multiplying 
on the right by the Pascal matrix P) we get to the matrix 



\ 



A 



( 1 
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1 














2 


2 


1 











4 


5 


3 


1 








9 


12 


9 


4 


1 





21 


30 


25 


14 


5 


1 



V 



/ 



It can be proved that A is the unique admissible matrix associated with 
Motzkin numbers. However, it will be an immediate consequence of the 
results of section El 



5.3 Bell numbers 

The most popular succession rule giving rise to the Bell numbers is the 
following: 

n -.\ (i)-(2) 

Such a rule describes the usual construction of set partitions, and its 
ECO matrix is the following: 



/ 1 

















• 


" ^ 





1 














• 







1 


1 











• 







1 


3 


1 








• 







1 


7 


6 


1 





• 







1 


15 


25 


10 


1 


• 







1 


31 


89 


65 


15 


1 • 




V : 
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The matrix F is the well-known matrix of the Stirling numbers of the 
second kind. Applying the fundamental change of basis leads to the matrix 

••• \ 
••• 
••• 

... • 

1 ••• 

\; ; \ --.J 

It is immediately seen that A is not an admissible matrix. Nevertheless, 
the elements of the main diagonal are all equal to 1, so A may be an Aigner 
matrix. Indeed, it can be shown that A is the Aigner matrix of type (er, r), 
where a = (£;)fc£N* an d r = (k — l)k>2 |A3j . This is our first example of an 
Aigner matrix which is not admissible and is linked to an ECO matrix by 
the fundamental change of basis. 



/ 1 











1 


1 








2 


3 


1 





5 


10 


6 


1 



15 37 31 10 



5.4 Factorial numbers 



The case of factorial numbers, which is extremely simple from the point 
of view of succession rules, turns out to be rather curious when the funda- 
mental change of basis is applied. Indeed, the trivial rule 



: 



(k) ~» (k + l) k 

for the factorial numbers leads to the diagonal ECO matrix 



/ 1 

10 

2 

6 

24 



V = 



\ 



/ 

where, clearly, / n & = (n — l)!<5 n fe {5 is the usual Kronecker delta). Thus we 
have P n (x) = (n - l)\x n = YJk=i( n ~ 1 ) ! (fc-i)^( a; )' whence 



A 



( 1 

110 

2 4 2 

6 18 18 6 

24 96 144 96 24 

V- ; : • ; 
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It is clear that A cannot be an Aigner matrix, since the elements on 
the main diagonal are not equal to 1. However, if we consider the scalar 
multiplication of the (n + l)-st and (m + l)-st rows of A, supposing that 
n < m, we get: 



n+l 

^ a (n+l)k a (m+l)k 
k=l 



= n\m\ 

k=0 



n\ I m 
k)\k 



n + m 
n 



l (n+m+l)l- 



(Equality W i s an application of Vandermonde's convolution). Thus A 
possesses a typical property of admissible matrices, without being neither 
admissible nor Aigner. 

5.5 Involutions 

Involutions are considered (from a succession rule point of view) in |FPj . 
They are generated by the following succession rule: 



Q : 



II) 



(k) (k — l) fc - 1 (A; H- 1) 
giving rise to the ECO matrix: 

\ 



/ 1 

















1 











1 





1 











3 





1 





3 





6 





1 



\ : : : : : •. / 
In this case, the fundamental change of basis leads to the matrix: 

\ 



/ 1 














1 


1 











2 


2 


1 








4 


6 


3 


1 





10 


16 


12 


4 


1 



V : 



/ 



which is the Aigner matrix of type (a, r), for a = (l)fcgN* and r = (k — l)k>2 
|A3j . This case has some analogies with that of Bell numbers (for example, 
A is Aigner but not admissible). 
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6 Factorial succession rules 

Referring to |B et a/. I IFF] , we recall the definition of a factorial succes- 
sion rule and a factorial rule operator. 

A factorial succession rule is a rule of the form: 



for k > ro > 1. A factorial rule operator is the rule operator of a factorial 
rule. In |B et al.\ it is shown that factorial rules have an algebraic generating 
function. In |FPj it is stated the following result, concerning the form of a 
factorial rule operator. 

Proposition 6.1 ()FP)[ ) A rule operator L is factorial if and only if L = 
p(x, x _1 , T), where p(a, b, c) is a polynomial of degree 1 in c having the form: 

p(a, b, c) = uo(a) + v (b) + ui(a)c, 

and T is the factorial derivative operator, as usual. 

In the present section we give a complete answer to the first of the 
problems stated in section 0] for the class of ECO matrices arising from fac- 
torial rules. In order to accomplish our result, we need to slightly generalize 
Aigner's original setting. 

Consider an infinite lower triangular matrix A = (a n k) n ,keN* > with an = 
1, and denote by L the linear operator associated with its rows, as in 0. We 
say that A is a generalized Aigner matrix when there exist three nonnegative 
integer sequences (r n ) ne N> («n)neN> (tn)neN such that, for every n € N: 



The following fact follows immediately from the above definition. 

Proposition 6.2 If A is a generalized Aigner matrix, then its entries obey 
the following recursion: 




(k) ~* (r )(r + 1) • • • (r + k - m - l)(k + d x ) • ■ • (k + d m ) 



L{p n (x)) = t n p n -l{x) + S n p n (x) + r n p n+ i(x). 



(10) 



I a (n+l)k — r fc-l a n(fc-l) + SkO-nk + tk+l a n(k+l) 

Proof. Since L{(5 n {x)) = (3 n+ \(x), using (jTU]l we have: 




(11) 



n+l 



n 




k=l 



k=l 



n 



^2 a n k (tkPk-i(x) + s k p k (x) + r k p k+1 (x)) 



k=l 
n+l 



^ (r fc _ia„( fc _i) + s k a nk + tjfe+ia n (fe+i)) Pk{x) 

k=l 
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whence the thesis follows. ■ 

In particular, it is easily seen that a generalized Aigner matrix is not 
forced to have the elements on the main diagonal equal to 1. 

Question: does the converse of proposition WM hold? 

Now we can focus on the case of factorial rules. Suppose that L is a 
factorial rule operator of the form 

L = a(x) + b^ 1 ) + c(x)T, (12) 

where a(x) = J2j t cikX k , b(x) = Ylk>i^kX k , c(x) = Ylk c k xk are fixed poly- 
nomials. We can immediately find a sufficient condition for L to induce a 
generalized Aigner matrix. 

Proposition 6.3 If b{x~ l ) = 0, deg a(x) < 1 and deg c(x) < 2, then A is 
a generalized Aigner matrix. 

Proof. We have to show that L acts as in (jlOjl . Indeed, a simple compu- 
tation shows that 

L(p n (x)) = ((a + aix) + (co + cix + c 2 x 2 )T)(p n {x)) 

= (co + ci + c 2 )p n -i(x) + (a + ai + ci + 2c 2 )p„(x) + (ai + c 2 )i?n 

which is enough to conclude thanks to the previous proposition. ■ 



Examples. 

i) The succession rules for Catalan and Motzkin numbers described above 
are associated with rule operators for which, respectively, a(x) = 
b[x) = 0, c(x) = x 2 and a(x) = x — 1, b(x) = 0, c(x) = x. 

ii) Consider the following succession rule, inducing Schroder numbers: 
(1) 

(1) - (2) . (13) 

(k)^(3)(4)(5)---(k-l)(k)(k + l) 2 



In this case, the rule operator L has the form: 

L = x - x 2 + x 3 T, 

so it does not satisfy the hypotheses of the above proposition. Never- 
theless, this rule is related to a generalized Aigner matrix, as we will 
see in the next pages. 
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Theorem 6.1 A factorial rule operator L as in is associated with a 
generalized Aigner matrix if and only if the following conditions hold: 

i) EkO a k + EkL k +1 )ck = 0, Vh>2; 
*) Ek^r-^r^h = 0, Vfc < n - 1. 
Proof. The rule operator L acts on the Aigner basis as follows: 

L(p n {x)) = ^ a k X k p n {x) + ^ b k X~ k p n {x) + ^ C k X k p n _i (x) 
k k>\ k 

k 



= Z^M Yl (l)Pn+h(x) J 
k \h=0 V J / 



k>l \h=0 

+ S C M Yl ( k h )pn-l+h{x) J • 
fc \h=0 V 7 / 

Now, L must satisfy condition ()1U|) . which means that, in the above 
expansion, all the coefficients of the polynomials Pk(x), for k ^ {n — 1, n, n + 
1}, must be zero. This translates into conditions i) and u) above, so the 
proof is complete. ■ 



Example. The previous example related to Schroder numbers can now 
be reconsidered. It is clear that condition ii) is trivially verified, whereas the 
only interesting case of condition i) occurs when h = 2, and we have: 



-1) + 



1 



-1 + 1 = 0. 



an 



So rule ()13|) is associated with a generalized Aigner matrix A. By 
explicit computation for the operator L, we find for A the following expres 



sion: 



/ 1 
1 

2 

6 

22 

V i 



o 
l 

4 
16 

68 





2 

14 

78 






4 

40 










\ 



where the entries obey the following recursion: 

- o„i + a n 2, 

- a n i + 3a„2 + a„3, 
= 2 a n(k-l) + 3anfc + a n (fc+i)> 




for k > 2 
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At this stage, it is worth noting that in [A3] an Aigner matrix for 
Schroder numbers is taken into consideration, precisely: 



A 



( 1 

2 1 

6 5 10 

22 23 8 1 

90 107 49 11 1 

\ ; ; ; ; : 



\ 



/ 



If we apply the inverse of the fundamental change of basis (that is, we 
switch from p n (x) to x n ), we find the following matrix: 

\ 



/ 1 














1 


1 











2 


3 


1 








6 


5 


10 


1 





22 


38 


22 


7 


1 



V : 



/ 



Strictly speaking, F is not an ECO matrix, as we can immediately no- 
tice. However, if we suppose that, in F, column k represents the distribution 
of the label (2k), we manage to find an ECO-interpretation. From a linear 
operator point of view, we can consider the operator L associated with the 
rows of F: if we replace the variable x with x 2 , we in fact obtain a rule oper- 
ator L, which corresponds to a well-known ECO-interpretation of Schroder 
numbers |BDT/PP2j . The succession rule related to L is the following: 



(2) 
(2*0 



(2)(4) 2 (6) 2 ---(2A;) 2 (2£; + 2) 



7 Differential succession rules 

We call differential succession rule each rule such that in the production 
of every label (k) at least one label greater or equal than (k — 1) has an 
exponent linearly depending on (k). Equivalently, a differential rule operator 
is a rule operator which can be expressed in the form 

L = p(x, x _1 , D) = o(x) + 6(x _1 ) + c(x)D. 

Observe that, for reasons of consistency (in a succession rule a node la- 
belled (k) must produce exactly k sons), in the above formula we necessarily 
have c(x) = x*, for some t G N, so that a differential rule operator has the 
following, general expansion: 

L = p(x, x~\ D) = o(x) + 6(x- x ) + x*D. (14) 
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The rules given in section E] for Bell numbers, factorials and the number 
of involutions of a set are examples of differential succession rules, whose 
associated differential rule operators are, respectively, x — 1 + xD, x 2 D and 
x — x _1 + D. In |FPj it is shown that the sequences determined by a differen- 
tial succession rule possess a transcendental (ordinary) generating function. 
Using arguments which are essentially analogue to those employed for fac- 
torial succession rules, we can prove the following results, concerning the 
relationship with Aigner's theory. 

Theorem 7.1 A differential rule operator L as in is associated with a 
generalized Aigner matrix if and only if the following conditions hold: 

1. if t ^ 0, the conditions are: 

i) Z k (t)ak + n( h t +1 )-{i-\)=0, Vh>2; 
H) £ fe >i(-l) n - fc -T^7% = °. Vh<n- 1; 

2. if t = 0, the conditions are: 

^E fe ffl«fc = 0, Vh>2; 

H) E^iC-ir-^C^K - = ' Vh<n- 1; 

Corollary 7.1 If b{x~ l ) = ; deg a{x) < 1 and t G {1,2}, then L is asso- 
ciated with a generalized Aigner matrix. 

Examples. The cases of Bell numbers and factorial numbers can be easily 
tackled using the corollary (see above for the rule operators involved). As 
far as involutions are concerned, we have to consider the rule operator x — 
x -1 + D. Applying theorem 17. II in the case t = 0, we have that condition i) 
is trivially satisfied, whereas condition ii) becomes: 

~ J ~ h \ (_i) _ (-I)"-* = {-\) n ~ h - (-l) n - h = 0, 
which is enough to conclude. 



8 Conclusions and further work 

The present work is intended to be only the first step towards a more 
detailed investigation of the relationship between the ECO method and the 
theory of Catalan-like numbers. 

One of the first things to be done in the next future is to provide a more 
complete gallery of examples to illustrate the soundness of our approach. 
In particular, it would be nice to find new applications of the ECO method 
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starting from known ones of Aigner's theory, and vice versa, which is what 
we hope to do in a forthcoming publication. 

Another line of research is provided by the last section of |A3j . where 
Aigner introduces the basics of what he calls "ballot enumeration". The 
analogies with the techniques of the ECO method are evident and, in fact, 
the combinatorial model proposed by Aigner is a particular instance of an 
ECO construction. 
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